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PROCEDE D ' ESTIMATION D'UN RAYONNEMENT DIFFUSE, 
NOTAMMENT AFIN DE CORRI6ER DES ME SURE S EN RADIOGRAPH-EE 

DESCRIPTION /^JX^ 
5 Le sujet de cette invention est un procede 

d' estimation d'un rayonnement diffuse, dont 
1 ' application principalement envisagee est la 
correction de radiograph ics . ^v***> ^ A^^y-t^ ^ **^***L*J£ 
L ' utilisati on . .f^'un _r ayonnement oon i ^ue , ^^^1 
10 tH?e« — fxag^^emre^ e^==^ad±ogra^*e , presente l'mconvenient ^ „ 

de produire un rayonnement diffuse important a travers < **^% t J 

>tA^-£^ r 1 'objet examine. En d'autres termes , chacun des 

/ ^ V, ^ ' 

- ^ i detecteurs situes derriere 1' objet recoit non seulement 
X f y-^r 

^^^sre- un rayonnement primaire, provenant directement de la 
^uJh 15 source par un trajet rectiligne et ayant traverse une 

^s^^Jz region bien definie de 1' objet, mais un rayonnement 

C C ^J^^j^ ^ diffuse de provenance indeterminee qui affecte la 
mesure et qu'il serait done souhaitable de corriger. 



Plusieurs procedes sont deja pratiques. 

20 C'est ainsi que le rayonnement primaire peut dtre 
mesure seul si une collimation stricte des detecteurs 
et de la source est faite afin d' intercepter le 
rayonnement diffuse , mais ce procede necessite en 
pratique un balayage du faisceau qui est lent a 

25 accomplir, et pendant lequel on doit s ' ac commode r de 
mouvements du patient si on examine des etres vivants. 

On a aussi eu l'idee contraire de ne 
mesurer que le rayonnement diffuse. On dispose pour 
cela un reseau discontinu d' absorbeurs , comme des 

30 billes de plomb, entre l'objet et les detecteurs, pour 
arreter localement le rayonnement primaire, de sorte 
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PROCEDE D' ESTIMATION D'UN RAYONNEMENT DIFFUSE, 
NOTAMMENT AFIN DE CORRIGER DES ME SURE S EN TOMOGRAPHIE 

OU OSTEDENSITOMETRIE 

5 DESCRIPTION 

Le sujet de cette invention est un procede 
d' estimation d'un rayonnement diffuse, dont 
1' application principalement envisagee est la 
correction de mesures en tomographie ou 

10 osteodensitometrie . 

L' utilisation d'un rayonnement 

d 1 irradiation, presente 1 ' inconvenient de produire un 
rayonnement diffuse important a t ravers l'objet examine 
surtout dans le cas tres frequent d ! un rayonnement 

15 divergent (conique ou en eventail) . En d'autres termes, 
chacun des detecteurs situes derriere l'objet re<?oit 
non seulement un rayonnement primaire, provenant 
directement de la source par un trajet rectiligne et 
ayant traverse une region bien definie de l'objet, mais 

2 0 un rayonnement diffuse de provenance indeterminee qui 
affecte la mesure et qu'il serait done souhaitable de 
corriger . 

Plusieurs procedes sont deja pratiques. 
C'est ainsi que le rayonnement primaire peut etre 

25 mesure seul si une collimation stricte des detecteurs 
et de la source est f aite af in d' intercepter le 
rayonnement diffuse, mais ce procede necessite en 
pratique un balayage du faisceau qui est lent a 
accomplir, et pendant lequel on doit s'accommoder de 

30 mouvements du patient si on examine des etres vivants. 



2 

que les detecteurs situes derriere ces absorbeurs ne 
mesurent que le rayonnement diffuse. Ce procede appele 
« beam stop » donne done des tables ou nappes 
bidimensionnelles de valeur de rayonnement diffuse , 
5 qu'on complete par interpolation entre les detecteurs 
places derriere les absorbeurs. Le rayonnement diffuse 
ainsi estime est soustrait du rayonnement total mesure 
separement. Ce procede est precis mais a 1 ' inconvenient 
qu'il impose deux irradiations de I'objet et done un 

10 doublement de la dose de rayons qu'il regoit. Un 
dernier exemple de methode de correction du rayonnement 
diffuse par des moyens materiels comporte l'emploi de 
grilles anti-dif f usantes , mais leur efficacite n'est 
que partielle ; elle est insuffisante pour un faisceau 

15 conique, ou le rayonnement diffuse peut etre plusieurs 
fois superieur au rayonnement primaire. 

Enfin, il existe un certain nombre de 
methodes numeriques pour estimer le rayonnement 
diffuse, a partir de convolutions ou de deconvolutions 

20 des mesures par exemple ; on pourrait aussi citer le 
brevet frangais 2 759 800 pour un procede numerique 
different , analytique. Elles sont en general d'emploi 
delicat car elles dependent de parametres choisis par 
1 'utilisateur (noyaux de convolution par exemple) qui 

25 ne donnent de bons resultats que dans des situations 
favorables, comme des petites zones ou le rayonnement 
diffuse est faible, ou des objets au contenu 
relativement homogene. Il n f existe aucun procede simple 
qui permette par exemple de corriger le rayonnement 

30 diffuse a travers le thorax ou d'autres grandes zones 
anatomiques, dont l'examen est frequent mais qui sont 

B 13^7^.3 JCI 
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On a aussi eu l'idee contraire de ne 
mesurer que le rayonnement diffuse. On dispose pour 
cela un reseau discontinu d' absorbeurs , comme des 
billes de plomb, entre l'objet et les detecteurs, pour 
5 arreter localement le rayonnement primaire, de sorte 
que les detecteurs situes derriere ces absorbeurs ne 
mesurent que le rayonnement diffuse. Ce procede appele 
« beam stop » donne done des tables ou nappes 
bidimensiormelles de valeur de rayonnement diffuse, 

10 qu'on complete par interpolation entre les detecteurs 
places derriere les absorbeurs. Le rayonnement diffuse 
ainsi estime est soustrait du rayonnement total mesure 
separement. Ce procede est precis mais a 1 ' inconvenient 
qu'il impose deux irradiations de l'objet et done un 

15 doublement de la dose de rayons qu'il regoit. Un 
dernier exemple de methode de correction du rayonnement 
diffuse par des moyens materiels comporte I'emploi de 
grilles anti-dif f usantes, mais leur efficacite n'est 
que partielle ; elle est insuffisante pour un faisceau 

2 0 conique, ou le rayonnement diffuse peut etre plusieurs 

fois superieur au rayonnement primaire. 

Enfin, il existe un certain nombre de 
methodes numeriques pour estimer le rayonnement 
diffuse, a partir de convolutions ou de deconvolut ions 
25 des mesures par exemple ; on pourrait aussi citer le 
brevet frangais 2 759 8 00 pour un procede numerique 
different, analytique. Elles sont en general d'emploi 
delicat car elles dependent de parametres choisis par 
l'utilisateur (noyaux de convolution par exemple) qui 

3 0 ne donnent de bons resultats que dans des situations 

favorables, comme des petites zones ou le rayonnement 
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defavorables pour corriger le rayonnement diffuse en 
raison de leur volume meme et de 1 ' heterogeneite due a 
la presence d'une structure d'os complexe et dont la 
capacite d' attenuation du rayonnement est tres 
5 differente de celle des tissus mous. 

Mentionnons enfin le brevet americain 
6 018 565 pour 1' expose d'une methode mixte, a « beam 
stop » et convolution. 

Un objet essentiel de 1' invention est de 
10 proposer un procede d' estimation et de correction de 
rayonnement diffuse qui puisse convenir pour des 
situations^ dif f iciles de radiographic ^ Q^r^A^^ ^ ^*>^^^ 
o( 9* * n ' ^ Le^proc6cfe ^confer ir^l*^ ' inv^^orTest / sou^^ 

sa forme la plus generale, un procede d'^s^^^i=o i\ A'xp& fi 
15 * i a^ u nn emea^— di^-fu-s- e prbvefm^t-d ' un rayonnement iftitxal 
"l^yarvt^t raver se un objet en subissant une attenuation 
laissant passer un rayonnement total de mesure, 
caracterise par : 

^ " I> une prise d'une table de mesures d'un 

2 0 rayonnement diffuse, obtenue en faisant passer le 
rayonnement initial par un simulacre de 1' objet, 
oJ<r ip**v <JLvt>*-e^^- — ■ Z s un calcul de coefficients de 

^S^^^transposition entre le simulacre et 1' objet , d'apres le 
rayonnement initial, le rayonnement total de mesure a 
25 travers 1' objet et un rayonnement total de mesure a 
travers le simulacre, 

- et une ponderation de la table de mesures 
avec les coefficients de transposition* 

Avantageusement , le simulacre sera un bloc 
30 d'epaisseur constante et en une matiere homogene, ayant 
une attenuation semblable a une matiere de base de 
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diffuse est faible, ou des objets au contenu 
relativement homogene . II n'existe aucun precede simple 
qui permette par exemple de corriger le rayonnement 
diffuse a travers le thorax ou d'autres grandes zones 
5 anatomiques, dont 1 ' examen est frequent mais qui sont 
defavorables pour corriger le rayonnement diffuse en 
raison de leur volume meme et de 1 ' heterogeneite due a 
la presence d'une structure d'os complexe et dont la 
capacite d' attenuation du rayonnement est tres 

10 differente de celle des tissus mous. 

Mentionnons enfin le brevet americain 
6 018 565 pour 1 ' expose d'une methode mixte, a « beam 
stop » et convolution. 

Un objet essentiel de 1' invention est de 

15 proposer un proc^de d' estimation et de correction de 
rayonnement diffuse qui puisse convenir pour des 
situations difficiles de controle non destructif 
d' objets inertes ou animes, ou tous procedes de 
reconstruction d f image. 

2 0 Le procede conforme a 1' invention est, sous 

sa forme la plus generale, un procede d'imagerie d f un 
objet par acquisitions multiples, comprenant une 
estimation d ! un rayonnement diffuse provenant d'un 
rayonnement initial ayant traverse un objet en 
25 subissant une attenuation laissant passer un 
rayonnement total de mesure, caracterise par : 

- pour au mo ins une acquisition, une prise 
d'une table de mesures d'un rayonnement diffuse, 
obtenue en faiisant passer le rayonnement initial par un 

3 0 simulacre de 1' objet, 
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l'objet ; en general la prise de table de mesure sera 
une selection dans une s<§rie de tables de raesures de 
rayonnement diffuse , obtenues auparavant en faisant 
successivement passer le rayonnement initial a travers 
5 une serie respective de simulacres de l'objet, 
d'epaisseurs differentes mais constante ; et la 
selection sera faite par comparaison d'une valeur du 
rayonnement total de mesure a travers l'objet et d'une 
valeur du rayonnement total de mesure a travers les 
10 simulacres. 

Les coefficients de ponderation sont 
g§neralement des rapports de valeurs d'une meme 
fonctionnelle calculee pour l'objet et pour le 
simulacre. La fonctionnelle utilisee peut etre egale au 

15 produit du rayonnement total de mesure par le 
logarithme du rapport de rayonnement total de mesure et 
du rayonnement initial. 

L' invention sera maintenant decrite en 
reference aux figures , parmi lesquelles : 

20 - la figure 1 est une vue generale d'une 

acquisition des mesures ; 

- la figure 2 est une vue d'une acquisition 
de calibration ; 

- et la figure 3 illustre les etapes du 

25 procede. 

Reportons-nous d'abord a la figure 1, ou un 
tube 1 de rayons X emet un faisceau 2 conique vers un 
objet 3 a examiner (ici un patient etendu sur une table 
4) puis, a travers lui, vers un reseau 5 plan de 
30 detecteurs 6 disposes en matrice. Les detecteurs 6 sont 
relies a un appareil d ' acquisition 7 et mesurent un 
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et pour chacune des acquisitions, un 
calcul de coefficients de transposition entre le 
simulacre et l'objet, d'apres le rayonnement initial, 
le rayonnement total de mesure a travers l'objet et un 
rayonnement total de mesure a travers le simulacre, 

- et une ponderation de la table de mesures 
avec les coefficients de transposition. 

Avantageusement , le simulacre sera un bloc 
d'epaisseur constante et en une matiere homogene , ayant 
une attenuation semblable a une matiere de base de 
l'objet ; en general la prise de table de mesure sera 
une selection dans une serie de tables de mesures de 
rayonnement diffuse, obtenues auparavant en faisant 
success ivement passer le rayonnement initial a travers 
une serie respective de simulacres de l'objet, 
d'epaisseurs differentes mais constante ; et la 
selection sera faite par comparaison d'une valeur du 
rayonnement total de mesure a travers l'objet et d'une 
valeur du rayonnement total de mesure a travers les 
simulacres . 

Les coefficients de ponderation sont 
generalement des rapports de valeurs d'une meme 
f onctionnelle calculee pour l'objet et pour le 
simulacre. La f onctionnelle utilisee peut §tre egale au 
produit du rayonnement total de mesure par le 
logarithme du rapport de rayonnement total de mesure et 
du rayonnement initial. 

L' invention sera maintenant decrite en 
reference aux figures, parmi lesquelles : 

- la figure 1 est une vue generale d'une 
acquisition des mesures ; 
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rayonneinent diffuse qui se superpose au rayonnement 
primaire, seul convenable pour La r^rilograpfrie^ 



L' estimation du rayonnement diffuse^ a 
tr avers le patient >i/conBxste tout d'abord a obtenir 
5 des tables bidimensionnelles ou nappes de rayonnement 
diffuse obtenues dans des circonstances comparables. 
Pour cela, on effectue des irradiations d'etalonnage a 
travers des simulacres 8 de l'objet 3^ a examiner, 
conform<§ment a la figure 2 : les conditions 
,10 d' irradiation restent les memes, c'est-a-dire qu'on 
continue d'utiliser le tube 1, le faisceau 2, le reseau 
5 de detecteurs 6 et l'appareil d' acquisition 7, le 
simulacre 8 remplagant cependant le patient ; on a 
aussi ajoute une grille 9 de billes 10 de plomb entre 
15 le simulacre 8 et le reseau 5* 11 resulte de cette 
disposition que les rayons 11 passant par les billes 10 
sont completement absorbes et que les regions 12 du 
reseau 5 situees dans le prolongement de ces rayons 11 
ont des detecteurs 6 qui ne mesurent que le rayonnement 
20 diffuse a ces endroits. 11 suffit de relever ces 
valeurs mesurees et d'interpoler entre les regions 12 
pour estimer convenablement le rayonnement diffuse issu 
du simulacre 8 pour tous les detecteurs 6 du reseau 5 . 

Le simulacre 8 devrait etre semblable a 
25 l'objet afin que les rayonnements diffuses par eux 
fussent identiques. Une similitude parfaite n'est pas 
realisable , et c'est pourquoi on se contente d'un 
simulacre 8 ressemblant a l'objet 3 et dont la nappe 
associee de rayonnement diffuse sera corrigee 
30 ulterieurement pour evaluer celle de l'objet. En 
pratique, le simulacre 8 peut etre un bloc d'une 
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- la figure 2 est une vue d'une acquisition 
de calibration ; 

- et la figure 3 illustre les etapes du 

procede . 

5 Reportons-nous d'abord a la figure 1, ou un 

tube 1 de rayons X emet un faisceau 2 conique vers un 
objet 3 a examiner (ici un patient etendu sur une table 
4) puis, £ travers lui, vers un reseau 5 plan de 
detecteurs 6 disposes en matrice. Les detecteurs 6 sont 
10 relies a un appareil d' acquisition 7 et mesurent un 
rayonnement diffuse qui se superpose au rayonnement 
primaire, seul convenable pour 1 1 examen ou le controle 
de 1 1 ob j et . 

L' estimation du rayonnement diffuse a 

15 travers le patient (1 ! objet 3) consiste tout d'abord a 
obtenir des tables bidimensionnelles ou nappes de 
rayonnement diffuse obtenues dans des circonstances 
comparables . Pour cela, on effectue des irradiations 
d'etalonnage a travers des simulacres 8 de 1'objet 3 a 

20 examiner, conformement a la figure 2 : les conditions 
d' irradiation restent les memes, c'est-a-dire qu'on 
continue d'utiliser le tube 1, le faisceau 2, le reseau 
5 de detecteurs 6 et 1' appareil d' acquisition 7, le 
simulacre 8 remplacpant cependant le patient ; on a 

25 aussi ajoute une grille 9 de billes 10 de plomb entre 
le simulacre 8 et le reseau 5. II resulte de cette 
disposition que les rayons 11 passant par les billes 10 
sont completement absorbes et que les regions 12 du 
reseau 5 situees dans le prolongement de ces rayons 11 

30 ont des detecteurs 6 qui ne mesurent que le rayonnement 
diffuse a ces endroits. II suffit de relever ces 
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matiere homogene et qui pr^sente le meme coefficient 
d' attenuation que la matiere de base de l'objet 3 : 
dans le cas d f un corps humain, compose pour l'essentiel 
de tissu mou, on sait que le plexiglas 
5 (polymethacrylate) convient. 

Afin de permettre des mesures varies, on 
disposera en realite de plusieurs nappes de rayonnement 
diffuse , obtenues pour autant de simulacres 8, qui ne 
differeront que par leur epaisseur et done par la 
10 longueur du trajet parcouru par les rayons 11* Ces 
nappes seront enregistrees dans une base de donnees 
prealablement aux mesures utiles sur les objets 3 
lihl i iiiijwaphicr - pour prendre une nappe de rayonnement 
diffuse comparable a celle d'un objet 3, on 
15 select ionnera en pratique une des nappes de la base de 
donnees ou, mieux, une nappe qu'on aura obtenue par des 
calculs d' interpolation entre deux de ces nappes. Le 
critere de selection pourra etre defini au moyen d'un 
rayon particulier 13 aboutissant a une region 14 du 
20 reseau 5 et qui ne passera ni par les absorbeurs 10 de 
la figure 2, ni par des tissus osseux du patient (ou 
plus generalement des portions de l'objet 3 dont les 
proprietes d ' absorption sont differentes du materiau du 
simulacre 8) a la figure 1. Le rayonnement total , 
25 primaire et diffuse / regu par la region 14 apres avoir 
traverse chaque simulacre 8 servira d' index a la table 
de rayonnement diffuse correspondante , et la table 
selectionnee aura 1' index a une valeur identique au 
rayonnement total mesure a la region 14 a tr avers 
30 l'objet 3. Tout cela correspond au passage de l'etat El 
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valeurs mesurees et d' interpoler entre les regions 12 
pour estimer convenablement le rayonnement diffuse issu 
du simulacre 8 pour tous les detecteurs 6 du reseau 5. 

Le simulacre 8 devrait §tre semblable a 
5 l'objet afin que les rayonnements diffuses par eux 
fussent identiques. Une similitude parfaite n'est pas 
realisable, et c'est pourquoi on se contente d'un 
simulacre 8 ressemblant a l'objet 3 et dont la nappe 
associee de rayonnement diffuse sera corrigee 

10 ulterieurement pour evaluer celle de l'objet. En 
pratique, le simulacre 8 peut etre un bloc d'une 
matiere homogdne et qui presente le meme coefficient 
d' attenuation que la matiere de base de l'objet 3 : 
dans le cas d'un corps humain, compose pour 1'essentiel 

15 de tissu mou, on sait que le plexiglas 
(polymethacrylate) convient . 

Afin de permettre des mesures variees, on 
disposera en r<§alite de plusieurs nappes de rayonnement 
diffuse, obtenues pour autant de simulacres 8, qui ne 

20 differeront que par leur epaisseur et done par la 
longueur du trajet parcouru par les rayons 11. Ces 
nappes seront enregistrees dans une base de donnees 
prealablement aux mesures utiles sur les objets 3. pour 
prendre une nappe de rayonnement diffuse comparable a 

25 celle d'un objet 3, on selectionnera en pratique une 
des nappes de la base de donnees ou, mieux, une nappe 
qu'on aura obtenue par des calculs d' interpolation 
entre deux de ces nappes. Le critere de selection 
pourra etre defini au moyen d'un rayon particulier 13 

30 about issant a une region 14 du reseau 5 et qui ne 
passera ni par les absorbeurs 10 de la figure 2, ni par 
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a 1'etat E2 dans 1 'organ i gramme de la figure 3, qu'on 
commence a coiranenter. 

La suite du procede consiste 

essentiellement en la correction de la table du 
5 rayonnement diffuse^ ainsi selectionne^pour l'ajuster 
au mieux qu'on puisse esperer a la nappe de rayonnement 
reellement diffuse par l'objet 3. Pour cela, on se sert 
de toutes les informations disponibles, c'est-a-dire du 
rayonnement total regu par les detecteurs 6 au-dela de 
10 l'objet 3 comme du siraulacre 8 selectionne* Ce 
rayonnement total etant note <E>t, le rayonnement diffuse 
<£d, le rayonnement initial issu du tube 1 <E>o et le 
rayonnement primaire <£, la relation 4>t = <I> + <&d est 
respectee. 

15 On est alors aux etats E3 et E4 de 

1 'organigramme de la figure 3. Ensuite, on transforme 
les valeurs des rayonnements totaux <E>t mesurees pour 
l'objet 3 et le simulacre 8 selectionne en leur 
appliquant des f onctionnelles . Plus precisement, il est 

2 0 connu dans l'art que <5d est proportionnel a 
<I>log(<I>/<I>o) ; cette relation, qui est appelee la loi de 
Klein et Nishina, donne une allure generale du 
rayonnement diffuse, h defaut de son intensite* 

Le rayonnement initial <I>0 est connu ; le 

25 rayonnement primaire <E> ne l'est pas, mais on consent a 
appliquer cette relation de fagon approchee en le 
remplagant par le rayonnement total <I>t, c'est-a-dire 
que la f onctionnelle employee associe a chaque valeur 
mesuree du rayonnement total <E>t la valeur calculee 

30 <J>tlog(Ot/4>o) , supposee proche du rayonnement diffuse 
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des tissus osseux du patient (ou plus generalement des 
portions de l'objet 3 dont les proprietes d' absorption 
sont differentes du materiau du simulacre 8) a la 
figure 1. Le rayonnement total, primaire et diffuse, 
regu par la region 14 apres avoir traverse chaque 
simulacre 8 servira d' index a la table de rayonnement 
diffuse correspondante, et la table select ionnee aura 
1' index a une valeur identique au rayonnement total 
mesure a la region 14 a travers l'objet 3. Tout cela 
correspond au passage de 1'etat El a l'etat E2 dans 
1 ' organigramme de la figure 3, qu'on commence a 
comment er . 

La suite du procede consiste 

essentiellement en la correction de la table du 
rayonnement diffuse ainsi selectionnee pour l'ajuster 
au mieux qu'on puisse esperer a la nappe de rayonnement 
reellement diffuse par l'objet 3. Pour cela, on se sert 
de toutes les informations disponibles, c'est-a-dire du 
rayonnement total regu par les detecteurs 6 au-dela de 
l'objet 3 comme du simulacre 8 selectionne. Ce 
rayonnement total etant note Ot, le rayonnement diffuse 
<3>d, le rayonnement initial issu du tube 1 Oo et le 
rayonnement primaire <D, la relation Ot - O + <I>d est 
respectee . 

On est alors aux etats E3 et E4 de 
1 ' organigramme de la figure 3. Ensuite, on transforme 
les valeurs des rayonnements totaux 4>t mesurees pour 
l'objet 3 et le simulacre 8 selectionne en leur 
appliquant des f onct ionnelles . Plus precisement, il est 
connu dans 1 ' art que <£d est proportionnel a 
Olog(O/0>o) ; cette relation, qui est appelee la loi de 



Od a cet endroit ; on est parvenu aux etats E5 et E6 de 
1 ' organ i gramme . 

L'etape suivante consiste a faire, pour 
chacun des detecteurs 6, le rapport des valeurs donnees 
5 par la f onctionnelle pour l'objet 3 et le simulacre 8 
selectionne selon la formule 

K = <>t Wgt / Oofrbjet ^ Les coeffici ents de 
<l>t log(<I>t / $o^imulacre 

ponderation K ainsi obtenus serviront a deformer la 

nappe de rayonnement diffuse selectionnee a 1'etat E2 

10 afin d'estimer celle de l'objet 3. Les resultats 
constituent encore une table bidimensionnelle ou une 
matrice ayant des dimensions identiques a celle des 
tables de rayonnement puisqu'elle est associee au 
reseau 5 de detecteurs 6. II est done possible et 

15 avantageux d'effectuer un filtrage numerique spatial de 
cette matrice en appliquant un filtre passe-bas qui 
corrige les coefficients K en ne conservant que les 
frequences les plus basses de leur variation et de les 
rendre ainsi probablement plus conformes a la realite 

2 0 puisque le rayonnement diffuse varie assez lentement 
d'un point a un autre. 

Quand la table des coefficients de 
ponderation definitifs, notes K' , a ete obtenue (a 
l'etat E7), elle sert a ponderer la table de 

25 rayonnement diffuse selectionnee auparavant a l'etat 
E2, pour obtenir une table de rayonnement diffuse par 
l'objet 3 (etat E8, qui constitue 1' estimation 
recherchee) ; la formule appliquee est 

<I>d objet=K' <Z>d simulacre, Ces valeurs estimees <J>d 

30 objet pourront alors etre soustraites du rayonnement 
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Klein et Nishina, donne une allure generale du 
rayonnement diffuse, a defaut de son intensite. 

Le rayonnement initial est connu ; le 

rayonnement primaire O ne l'est pas, mais on consent a 
appliquer cette relation de fagon approchee en le 
remplagant par le rayonnement total Ot, c'est-a-dire 
que la f onctionnelle employee associe a chaque valeur 
mesuree du rayonnement total <I>t la valeur calculee 
<t>tlog (Ot/Oo) , supposee proche du rayonnement diffuse 
<&d a cet endroit ; on est parvenu aux etats E5 et E6 de 
1 ' organigramme . 

L'etape suivante consiste a faire, pour 

chacun des detecteurs 6, le rapport des valeurs donnees 

par la f onctionnelle pour l'objet 3 et le simulacre 8 

selectionne selon la formule 

v - Ot log(g>t / <E>o>>bjet ^ . . 

K " 4>t log(*t / Oo>imulacre ' LeS coefficients de 

ponderation K ainsi obtenus serviront a deformer la 
nappe de rayonnement diffuse selectionnee a 1'etat E2 
afin d'estimer celle de l'objet 3. Les resultats 
constituent encore une table bidimensionnelle ou une 
mat rice ayant des dimensions identiques a celle des 
tables de rayonnement puisqu'elle est associee au 
reseau 5 de detecteurs 6. II est done possible et 
avantageux d'effectuer un filtrage numerique spatial de 
cette matrice en appliquant un filtre passe-bas qui 
corrige les coefficients K en ne conservant que les 
frequences les plus basses de leur variation et de les 
rendre ainsi probablement plus conformes a la realite 
puisque le rayonnement diffuse varie assez lentement 
d'un point a un autre. 
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total Ot mesure par les detecteurs 6 pour estimer le 
rayonnement primaire O et obtenir une image 
radiographique plus precise de l'objet 3. ^ , 

Ce procede s' applique aux ^^Ljj^x^hxB^ a 
energie >v jrhR^T ^n simple ou multiple ; dans le 
second cas, il est repete separ^ment pour chacune des 
energies employees . 

La fonctionnelle proposee ici n'est pas la 
seule qu'on puisse employer, et la fonctionnelle plus 
simple <I>d=k<I> (approchee ici encore en ^d^kOt), k etant 
une constante, pourrait aussi donner de bons resultats 
pour estimer <l>d. 
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Quand la table des coefficients de 
ponderation def initifs, notes K' , a ete obtenue (a 
l'etat E7) , elle sert a ponderer la table de 
rayonnement diffuse selectionnee auparavant a l'etat 
5 E2, pour obtenir une table de rayonnement diffuse par 
l'objet 3 (etat E8 , qui constitue l'estimation 
recherchee) ; la formule appliquee est 

<3>d objet=K' <Dd simulacre. Ces valeurs estimees <l>d 
objet pourront alors etre soustraites du rayonnement 

10 total <&t mesure par les detecteurs 6 pour estimer le 
rayonnement primaire et obtenir une image 

radiographique plus precise de 1' objet 3. 

Ce procede s' applique aux irradiations a 
energie simple ou multiple ; dans le second cas, il est 

15 repete separement pour chacune des energies employees. 

La f onctionnelle proposee ici n'est pas la 
seule qu'on puisse employer, et la f onctionnelle plus 
simple <I>d=k<I> (approchee ici encore en <£d=k<&t) , k etant 
une constante, pourrait aussi dormer de bons resultats 

20 pour estimer 3>d. 

Ainsi qu'on l'a mentionne, ce procede est 
particulierement interessant en tomographie ou 
osteodensitometrie, ou une image en profondeur de 
l 1 objet est reconstruite a partir d'une multitude 

2 5 d* acquisitions (irradiations) prises autour de 1' objet 
par un reseau mobile de detecteurs, puis par une 
combinaison numerique de ces acquisitions (menee par 
des techniques qui ne font pas 1' objet de 1* invention) . 
Une des limitations rencontrees en pratique est la dose 

30 excessive regue par l 1 objet ; or on a vu que le procede 
conforme a l f invention permettait de diminuer la dose 
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necessaire par rapport a d'autres, surtout ceux qui 
imposent une double irradiation pour soustraire le 
rayonnement diffuse du rayonnement total. On pressent 
que cet avantage sera sensible dans les procedes a 
5 reseau lineaire de detecteurs (ou a lignes superposees 
dans un reseau bidimensionnel ou on obtient une ou 
plusieurs coupes (images bidimensionnelles) a travers 
1 1 objet . 

Quand plusieurs acquisitions k sont 
10 entreprises, le mode operatoire resume par la figure 3 
peut etre applique a chacune d'elles, une estimation 
independante du rayonnement diffuse etant faite a 
chaque fois. II existe pourtant des situations ou le 
rayonnement diffuse pourra etre juge invariable, 
15 notamment pour des objets a symetrie de revolution, 
tres frequente en procedes de controle non destructif. 
Plusieurs des etapes de la figure 3 deviendront alors 
inutiles : on pourra ainsi se contenter de la selection 
d'un seul simulacre 8 pour toutes les acquisitions k ; 
2 0 event uellement , seules les etapes E3, E5, E7 et E8 
impliquant 1' objet devront etre r^petees pour chaque 
acquisition. Le choix d'une simplification du procede 
par 1 ! omission de certaines calibrations pourra etre a 
la discretion de l'operateur. 




25 
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«i — ^i^^k: ~^ju^x, 

REVEND I CATIONS^ 

1. Procede $$' estimation d'un rayonnement 
diffuse provenant d'un rayonnement initial ayant 
traverse un objet (3) en subissant une attenuation 
5 laissant passer un rayonnement total de mesure, 

A**. ^ caracterise par : 

^ — IT" une prise d'une table de mesures d'un 

^ ^^^w, rayonnement diffuse, obtenue en faisant passer le 
rayonnement initial par un simulacre (8) de l'objet, 
10 ^ un calcul de coefficients (K' ) de 



.*)r ^ transposition entre le simulacre et 1' objet f d'apres le 

^^A^ rayonnement initial (<E>o) , le rayonnement total de 

mesure a travers 1' objet (<E>t objet) et un rayonnement 
total de mesure a travers le simulacre (Ot simulacre), 
15 - et une ponderation de la table de mesures 

avec les coefficients de transposition. ^ 

2, Procede d' estimation — d-Ua£- — ^^ayenrrem^nt 
dif^^se selon la revendication 1, caracterise en ce que 
le simulacre (8) est un bloc d'epaisseur constante et 

2 0 en une matiere homogene, ayant une attenuation 
semblable a une matiere de base de 1' objet. 

3. Procede d' e-s^ima -cH-tm — rayonnom cnt - 
dif fuse selon la revendication 1, caracterise en ce que 
la prise de table de mesures est une selection dans une 

25 serie de tables de mesures de rayonnement diffuse, 
obtenues en faisant successivement passer le 
rayonnement initial a travers une serie respective de 
simulacres de 1' objet, qui sont des blocs d'epaisseurs 
differentes mais constante et en une matiere homogene, 

30 ayant une attenuation semblable a une matiere de base 
de 1' objet. 
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RBVEND I CAT I ONS 

1. Procede d 1 imager ie d'un objet par 
acquisitions multiples, comprenant une estimation d'un 
rayonnement diffuse provenant d'un rayonnement initial 
5 ayant traverse un objet (3) en subissant une 
attenuation laissant passer un rayonnement total de 
mesure, caracterise par : 

- pour au moins une des acquisitions, une 
prise d'une table de mesures d'un rayonnement diffuse, 

10 obtenue en faisant passer le rayonnement initial par un 

simulacre (8) de 1' objet, 

et pour chacune des acquisitions, un 

calcul de coefficients (K' ) de transposition entre le 

simulacre et 1' objet, d'apres le rayonnement initial 
15 (Oo) , le rayonnement total de mesure a travers 1' objet 

(<3>t objet) et un rayonnement total de mesure a travers 

le simulacre (Ot simulacre) , 

- et une ponderation de la table de mesures 
avec les coefficients de transposition. 

2 0 2. Proclde d' imagerie selon la 

revendication 1, caracterise en ce que le simulacre (8) 
est un bloc d'epaisseur constante et en une matiere 
homogerie, ayant une attenuation semblable a une matiere 
de base de 1' objet. 
25 3. Procede d' imagerie selon la 

revendication 1, caracterise en ce que la prise de 
table de mesures est une selection dans une serie de 
tables de mesures de rayonnement diffuse, obtenues en 
faisant success ivement passer le rayonnement initial a 

3 0 travers une serie respective de simulacres de 1' objet, 

qui sont des blocs d'epaisseurs differentes mais 
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4. Procede d ' ^fit i m** i on " n ray onn o m o f rb 

diffus e selon la revendication 3, caracterise en ce que 
la selection comprend une interpolation entre deux des 
tables de mesures. 



5 5. Procede d' as±±ma±j — d-M*a — ^aye^imeme^nir 

di££«-se selon l'une quelconque des revendications 3 et 
4, caracterise en ce que la selection est faite par 
comparaison d'une valeur du rayonnement total de mesure 
a travers l'objet et d'une valeur du rayonnement total 

10 de mesure a travers les simulacres. 

6. Procede d' Q E> tima - tion — d ' uffi^r a yomementr- 
diffu se selon la revendication 5, caracterise en ce que 
la comparaison est faite pour des rayons identiques 
(13) du rayonnement initial a travers l'objet et les 

15 simulacres, ne traversant que la matiere de base de 
l'objet. 

7 • Procede d ' ^sti-ma^ioff — d~^un- — r-^yoiinemejxt^^ 
^i#fnse selon l'une quelconque des revendications 1 a 
6, caracterise en ce que les coefficients de 

2 0 ponderation sont des rapports de f onctionnelle 
identiques calculees pour l'objet ^et pour le simulacre. 

8 . Procede d ' esL ima t i el^^ aa — -rayon ugmgjit _ 
rji ffngp selon la revendication 7, caracterise en ce que 
les fonctionnelles sont egales au produit du 

25 rayonnement total de mesure par le logarithme du 

rapport du rayonnement total de mesure et du 

rayonnement initial. 

9. Procede d ' es^timarti^m- — d-M**i — rayeR-nement 

^difiue-e selon l'une quelconque des revendications 1 a 
30 8, caracterise en ce qu'il comprend une etape de 

filtrage passe-bas des coefficients de transposition, 
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constante et en une matiere homogene, ayant une 
attenuation semblable a une matiere de base de l'objet. 

4. Procede d'imagerie selon la 
revendication 3, caracterise en ce que la selection 

5 comprend une interpolation entre deux des tables de 
mesures . 

5. Procede d'imagerie selon l'une 
quelconque des revendications 3 et 4, caracterise en ce 
que la selection est faite par comparaison d'une valeur 

10 du rayonneraent total de mesure a travers l'objet et 
d'une valeur du rayonnement total de mesure a travers 
les simulacres . 

6. Procede d'imagerie selon la 
revendication 5, caracterise en ce que la comparaison 

15 est faite pour des rayons identiques (13) du 
rayonnement initial a travers l'objet et les 
simulacres, ne traversant que la matiere de base de 
1 ' ob j et . 

7. Procede d'imagerie selon l'une 
20 quelconque des revendications 1 a 6 f caracterise en ce 

que les coefficients de ponderation sont des rapports 
de f onctionnelle identiques calculees pour l'objet et 
pour le simulacre. 

8. Procede d'imagerie selon la 
25 revendication 7, caracterise en ce que les 

f one tionne lies sont egales au produit du rayonnement 
total de mesure par le logarithme du rapport du 
rayonnement total de mesure et du rayonnement initial. 

9. Procede d'imagerie selon l'une 
30 quelconque des revendications 1 a 8, caracterise en ce 

qu'il comprend une etape de filtrage passe-bas des 



arranges en une table superposable a la table de 
mesures . 

10. Procede d* — a^adieg-^aplii-e comprenant une 
etape de correction de mesures de radiographie par une 
soustraction d'un rayonnement diffuse estime selon le 
procede selon l'une quelconque des revendications 
precedentes . 
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coefficients de transposition, arranges en une table 
superposable a la table de mesures . 

10. Procede comprenant une etape de 
correction de mesures de radiographie par une 
5 soustraction d'un rayormement diffuse estime selon le 
procede selon l'une quelconque des revendications 
precedentes . 

11) Application du procede selon l'une 
quelconque des revendications precedentes a la 

10 tomographie. 

12) Application du procede selon l'une 
quelconque des revendications 1 a 11 a 
1 ! osteodensitometrie . 

13) Application du procede selon l'une 
15 quelconque des revendications 1 a 11 au controle non 

destructif . 
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INPI 



I INSTITU1 
NATIONAL DE 
LA PROPMETE 
INDU5TRIEILC 



BREVET D'INVENTION 

CERTIFICAT D'UTILITE 

Code de la propriete intellectuelle - Livre VI 



N° 11235*02 



D&PARTEMENT DES BREVETS 

26 bis, rue de Saint P6tersbourg 
75800 Paris Cedex 08 

Telephone : 01 53 04 53 04 Telecopie : 01 42 93 59 30 



DESIGNATION D'INVENTEUR(S) Page N° .1. . / J. . 

(Si le demandeur n'est pas I'inventeur ou I'unique inventeur) 



Cet im prime est a remplir lisiblement a I'encre noire 



DB 113W/26G899 



Vos references pour ce dossier 

(facultatif 


B 14164.3/JCI DD2365 


N° D'ENREGISTREMENT NATIONAL 


02.10422 du 20.08.2002 



UTRE DE L'INVENTION (200 caracteres ou espaces maximum) 

PROCEDE D'ESTIMATION D'UN RAYONNEMENT DIFFUSE, NOTAMMENT AFIN DE CORRIGER 
DES MESURES EN TOMOGRAPHIE OU OSTEDENSITOMETRIE 



LE(S) DEMANDEUR(S) : 

COMMISSARIAT A L'ENERGIE ATOMIQUE 
3 1/33 rue de la Federation 
75752 PARIS 15eme 



DESIGNE(NT) EN TANT QU'INVENTEUR(S) : (Indiquez en haut a droite «Page N° 1/1» S'il y a plus de trois inventeurs, 
utilisez un formulaire identique et numerotez chaque page en indiquant !e nombre total de pages). 



Nom 


DINTEN 


Prenoms 


Jean-Marc 


Adresse 


Rue 


138, avenue des Freres Lumiere 


Code postal et ville 


69008 | LYON 


Societe d'appartenance (facultatif) 




Nom 


DARBOUX 


Prenoms 


Michel | 


Adresse 


Rue 


17 rue Aime Berey 


Code postal et ville 


38000 GRENOBLE 


Societe d'appartenance (facultatif) 




Nom 




Prenoms 




Adresse 


Rue 




Code postal et ville 




Societe d'appartenance (facultatif) 




DATE ET SIGNATURE(S) 
DU (DES) DEMANDEUR(S) 
OU DU MANDATAIRE 
(Nom et qualite du signataire) 

PARIS LE 1 1 Septembre 2002 J 
P. RICHARD C> Y_\ (jJ^C^^ 
422-5/002 (^^^^ 





La loi n°78-17 du 6 janvier 1978 relative a Tinformatique, aux fichiers et aux libertes s'applique aux reponses faites a ce formulaire. 
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